AoB PLANTS

Auxins are a family of plant hormones involved in various processes across plant tissues and species. The Nuclear Auxin Pathway (NAP) consists of interacting transcription factors (ARFs) and repressors (Aux/IAAs), which govern an individual cells response to changes in auxin concentration. These components are present in all land plants, and many species possess multiple copies of each signalling component. We present a general framework for ODE-based models of NAP submodules with the flexibility to model the promotion and repression of target genes by any combination of transcriptional regulators. We analyse published data and show that auxin treatment in Arabidopsis thaliana roots triggers a range of characteristically distinct temporal response profiles--for both target genes and the signalling components themselves. Using our modelling framework, we recapitulate aspects of this behaviour by presenting examples of real and theoretical NAP subnetworks, and by analysing the effect that these network dynamics have on auxin-mediated transcriptional responses. This work demonstrates the utility of our modelling framework as a general-purpose tool for understanding the function of certain protein-protein and protein-DNA interactions through their effects on the NAP. This exploration of the rich dynamics of more complex signalling pathways promises to advance our understanding of the NAP.

Soybean, a widely cultivated leguminous crop valued for its protein, amino acids, and oil, faces the challenge of maintaining protein levels, which have an inverse correlation with yield. Reducing leaf chlorophyll levels could increase seed protein levels without compromising yield; however, this is yet to be tested. Therefore, to understand the impacts of low chlorophyll mutations on soybean yield and seed composition, we screened and compared 25 low chlorophyll soybean mutants to their 11 dark green parents. PI548210 (Lincoln mutant) demonstrates a higher concentration of protein without affecting yield compared to its dark green parent PI548362 (Lincoln), suggesting it as a good candidate for further large-scale field trials. PI547555 (Y11/y11, Clark mutant) demonstrates a lower concentration of oil without impacting yield, alongside lower gross photosynthesis, but with chlorophyll levels in the pod and seed tissues that are comparable to its dark green parent PI548533 (Clark). These findings are consistent with the oil concentration of the soybean being influenced by pod and seed photosynthesis, which is correlated with pod height and row spacing. Chlorophyll levels in the leaf do not necessarily correlate with those in the pod and seed of low chlorophyll mutants, possibly due to substantially lower expression of chlorophyll synthesis genes in the pod and seed. SIGNIFICANCEO_LIPI548210 (Lincoln mutant), one of twenty-five low chlorophyll soybean mutants, demonstrates a higher concentration of soybean protein without affecting yield compared to its dark green parent (Figure 1 and Table 1). C_LIO_LIPI547555 (Y11/y11, Clark mutant), a low chlorophyll soybean mutant, demonstrates a reduced concentration of soybean oil without impacting yield, alongside lower gross photosynthesis in pod and seed tissues compared to its dark green parent (Figures 3 and Table 2). These findings suggest that the oil concentration of the soybean is influenced by pod and seed photosynthesis, which is in turn influenced by pod height and row spacing (Figure 2). C_LIO_LIChlorophyll levels in the leaf do not necessarily correlate with those in the pod and seed of low chlorophyll mutants, possibly due to substantially lower expression of chlorophyll synthesis genes in the pod and seed (Figure 5-6). C_LI O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=84 SRC="FIGDIR/small/744892v1_fig1.gif" ALT="Figure 1"> View larger version (55K): org.highwire.dtl.DTLVardef@4282dcorg.highwire.dtl.DTLVardef@9d565forg.highwire.dtl.DTLVardef@1918292org.highwire.dtl.DTLVardef@1359b1_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 1.C_FLOATNO Two low chlorophyll mutants are as healthy as their dark green parents. Lincoln and its low chlorophyll mutant, left; Clark and its low chlorophyll mutant, known as Y11/y11, right. It can be seen by eye that the plants have low chlorophyll (light green/yellow leaves) but a similar growth habit to their dark green parents. See Supplemental Figures 1-4 for contrast, where low chlorophyll mutants are stunted in growth compared to their dark green parents. C_FIG O_TBL View this table: org.highwire.dtl.DTLVardef@657ec9org.highwire.dtl.DTLVardef@166e75borg.highwire.dtl.DTLVardef@df23c7org.highwire.dtl.DTLVardef@1a60124org.highwire.dtl.DTLVardef@194ed96_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 1.C_FLOATNO O_TABLECAPTIONComparison of seed yield, weight, seed composition between low chlorophyll mutants and their dark green parents. ANOVA is used with linear mixed model (random effect = block, fixed effect = variety). Least squares mean is used to compare. For yield and seed composition, N=4 blocks. For leaf chlorophyll (SPAD), N=40. Yield is average yield per plant (g). n.s. = not significant. C_TABLECAPTION C_TBL O_FIG O_LINKSMALLFIG WIDTH=179 HEIGHT=200 SRC="FIGDIR/small/744892v1_fig3.gif" ALT="Figure 3"> View larger version (26K): org.highwire.dtl.DTLVardef@7a368aorg.highwire.dtl.DTLVardef@192b8f0org.highwire.dtl.DTLVardef@1abb738org.highwire.dtl.DTLVardef@89e978_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 3.C_FLOATNO Light response curve of low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533). Rates of net and gross photosynthesis of low chlorophyll (white) and dark green parents (black) pods under field conditions. Each dot represents a value (n=4) {+/-}SE. We assumed that the seeds greatly inhibited the transmittance of light through the pod and used photosynthetic photon flux density for a single-side. C_FIG O_TBL View this table: org.highwire.dtl.DTLVardef@3f0528org.highwire.dtl.DTLVardef@16ba712org.highwire.dtl.DTLVardef@a5ab2aorg.highwire.dtl.DTLVardef@889254org.highwire.dtl.DTLVardef@3efa4f_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 2.C_FLOATNO O_TABLECAPTIONPod photosynthetic parameters for low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533). Photosynthesis was measured 1 September through 15 September 2021 at the University of Illinois Energy Farm in Urbana, IL, USA. The statistical analysis was done using ANOVA with linear mixed model (alpha=0.05). N=4 {+/-} SEM for Clark and N=3 {+/-} SEM for Y11. C_TABLECAPTION C_TBL O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=130 SRC="FIGDIR/small/744892v1_fig2.gif" ALT="Figure 2"> View larger version (23K): org.highwire.dtl.DTLVardef@a36c26org.highwire.dtl.DTLVardef@1116c8forg.highwire.dtl.DTLVardef@ee5e61org.highwire.dtl.DTLVardef@1766712_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 2.C_FLOATNO Low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533) differ in concentration of seed oil, which interacts with height of pod and row spacing. The box plots show the median (central line), the lower and upper quartiles (box) and the minimum and maximum values (whiskers). The statistical analysis was done using ANOVA with linear mixed model (n=3 blocks, alpha=0.05). Least squares mean is used to compare. N.s., non- significant in the analysis. A. Concentration of oil in low chlorophyll mutant seeds from the upper canopy decreased by 4% compared to the dark green parent (18.2% vs 19%) while there was no difference between them in the seeds from the lower canopy (20.2% vs 20.6%). B. Schematic layout of 2013 field setting showing two different row spacings. C. Concentration of oil in low chlorophyll mutant decreased by 2% in 38cm spacing (21.4% vs 22%) while there was no difference in 19cm spacing (21.3% vs 21.7%) in 2013 field. C_FIG O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=162 SRC="FIGDIR/small/744892v1_fig5.gif" ALT="Figure 5"> View larger version (22K): org.highwire.dtl.DTLVardef@68e508org.highwire.dtl.DTLVardef@94a6ccorg.highwire.dtl.DTLVardef@152a187org.highwire.dtl.DTLVardef@1eae137_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 5C_FLOATNO (greenhouse). Correlation between the level of leaf chlorophyll (x-axis: SPAD reading) and the level of immature pod or seed chlorophyll (y-axis, mg/g DW). Line represents the linear regression model. R-squared is a coefficient of determination, the percentage of the response variable variation that is explained by the linear model. Pod is labeled by the fresh weight of seeds it contained. A. Level of chlorophyll of 25-100mg pod (n=18). B. Level of chlorophyll of 100-200mg pod (n=17) . C. Level of chlorophyll of 25-100mg seed (n=17). D. Level of chlorophyll of 100-200mg seed (n=20). C_FIG O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=180 SRC="FIGDIR/small/744892v1_fig6.gif" ALT="Figure 6"> View larger version (28K): org.highwire.dtl.DTLVardef@167fd88org.highwire.dtl.DTLVardef@361472org.highwire.dtl.DTLVardef@786325org.highwire.dtl.DTLVardef@1b53855_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 6.C_FLOATNO Levels of gene expression in chlorophyll synthesis pathway. A. CHL common pathway genes; Glutamyl-tRNA reductase (GluTR). Glutamate 1- semialdehyde aminotransferase (GSA-AT). ALA dehydratase (ALAD). Uroporphyrinogen III synthase (UROS). Uroporphyrinogen III decarboxylase (UROD). Protoporphyrinogen IX oxidase (PPO). B. Mg branch; Mg-chelatase (Mgch). Magnesium-protoporphyrin IX monomethyl ester cyclase (MPEC). Protochlorophyllide reductase (POR). 3,8-divinyl protochlorophyllide a 8-vinyl-reductase (4VCR). Heme pathway; Ferrochelatase (FECH). Heme oxygenase (HO). Phytochromobilin synthase (HY). Data come from Severin et al (2010). RPKM, reads per kilobase per million mapped reads. DAF, days after flowering. The source seed is experimental line A81-356022 which was generated by introgressing G. soja (PI468916) into G. max (A81-356022). C_FIG

A strategy to improve agricultural water productivity is to increase water use efficiency (WUE) at the level of plant transpiration through genetic selection. This requires detectable genetic variability in WUE and the ability to phenotype and select plants with higher WUE within a population. A proxy for phenotyping leaf level WUE by measuring carbon isotope signature ({delta}13Cleaf) has been supported by theory and data in C4 species. However, the functional relationship of {delta}13Cleaf and WUE in C4 species can be driven by genetics and environment. Therefore, a wide survey of existing natural variation is needed to quantify the heritability and identify various genetic factors that influence {delta}13Cleaf and WUE. In this study a genome-wide association panel was used to quantify the heritability of {delta}13Cleaf. We measured {delta}13Cleaf across a population of 360 genetically diverse lines of the C4 species Sorghum bicolor with single nucleotide polymorphic (SNP) markers determined from whole-genome resequencing. This analysis was conducted on two independent field environments where heritability of {delta}13Cleaf was evident and was driven by small genetic effects from loci that were consistently identified across environments. Candidate genes are presented that offer insights on future targets to manipulate and explore the functional relationship between {delta}13Cleaf and WUEi in C4 plants.

Cannabis is well known for its pungent, skunk-like aroma. Recent chemical studies have identified prenylated and C6 volatile sulfur compounds as contributors to its skunky and citrus-like aromas, but the pathways that produce these compounds remain unknown. This gap limits efforts to explain variation in sulfur-aroma traits and to selectively enhance or reduce those traits. To address this gap, we used the known chemistry of sulfur-containing volatiles in Cannabis and characterized sulfur and volatile biosynthetic pathways in other plant species to select candidate enzyme groups. Because the GMO cultivar is anecdotally associated with a pronounced sulfurous aroma, reference protein sequences and profile hidden Markov models were used to search its version 1 (v1) primary high-confidence protein set of 55,790 sequences. These searches recovered 975 unique proteins. Sequence screening retained 941 candidates across 20 reporting categories; 939 contained all expected domains, while the two candidates assigned to the methionine gamma-lyase (MGL)-nearest category had no category-specific expected-domain rule. The largest reporting category comprised 359 proteins containing a cytochrome P450 domain, recovered through a search motivated by cytochrome P450 family 74 (CYP74) enzymes involved in oxylipin and plant volatile formation. Thirteen of these proteins were also recovered by at least one full-length CYP74 reference search. Other large reporting categories included 218 sugar-transferase, 83 glutathione-transferase, and 61 alcohol dehydrogenase candidates. Comparison with the Cannabis Expression Atlas linked 168 candidates to 128 annotated genes through 100%-identity amino-acid matches spanning at least 80% of each GMO v1 candidate protein. Twenty-nine genes were tissue-specific, including 13 root-specific and 6 trichome-specific genes. These results define candidates for biochemical testing and direct searches for additional enzymes acting upstream and downstream in Cannabis sulfur-volatile pathways.

Blue Economy models position aquaculture as a key pathway to securing global food security. Species selected for aquaculture typically show rapid growth, high stress tolerance and fast biomass accumulation, but these same traits may increase their potential to become invasive when introduced beyond their native range. We investigated the invasion history and current status of the commercially important red seaweed Kappaphycus alvarezii in the Palk Bay-Gulf of Mannar region of India. This is one of the worlds largest cultivation hubs, a climatically vulnerable marine biodiversity hotspot, and one of the three regions to report invasion. We combined in-water surveys, interviews with wild seaweed collectors, and a review of published literature to reconstruct the history of invasion and assess current status. Invasion has declined substantially, with interviews indicating that the disappearance of invasive populations began around 2014. We discuss several non-mutually exclusive explanations for this decline, including climate change, loss of coral substrate, herbivory, and reduced vitality of the seaweed. Although the decline in invasion is encouraging for coral reefs, our findings raise questions about the ecological and socioeconomic consequences of introducing non-native aquaculture species under Blue Economy initiatives, particularly in ecologically sensitive regions vulnerable to climate change.

Plants competing belowground may produce extra roots, fewer roots, or no detectable change compared to plants growing alone. This inconsistency is often attributed to plants altering their root allocation in response to diverse cues, including neighbor detection and resource depletion by neighbors, but isolating these cues experimentally is challenging. Here, we hypothesize that water depletion alone can generate the range of root allocation strategies reported in the literature. We present this hypothesis as a water-explicit optimization model of root allocation that predicts a non-monotonic response. The model identified a critical depletion rate at which allocation shifted from increasing to decreasing with depletion. We tested this prediction using artificially rooted pots that imposed controlled water depletion while excluding living neighbors and their cues. A continuous artificial depletion gradient revealed the predicted hump-shaped pattern. These results reframe root overproliferation and underproliferation as positions along a single depletion-response curve.

Root hairs (RH) are excellent model systems for studying cell size and polarity since they elongate several hundred-fold their original size. Their tip growth is regulated by both intrinsic and environmental signals and is associated with the existence of a highly controlled cytoplasmic tip Ca{superscript 2} gradient, whose disruption impairs RH development. The molecular mechanisms underlying the Ca2+ homeostasis fine tuning and the Ca2+ organellar contributions to the cytoplasmic pool remain unclear. In the model plant Arabidopsis thaliana, many efflux routes are present, including those that employ Ca2+-pumps from the Autoinhibited Ca2+-ATPase (ACA) family. Here, we identified that the ER localized ACA2, and to a lower extent ACA7, are crucial ACAs required to control RH growth. By using genetically encoded Ca2+ biosensors we showed that Ca2+-dynamics are compromised in the aca2-2 mutant, having lower cytosolic Ca2+ concentration [Ca2+]cyt and growth rate, showing an altered homeostatic calcium setpoint compared to Col-0. Accordingly, the ACA2 mutation changed the dynamics of [Ca2+]cyt oscillations coupled to growth rate, inducing longer periods and more regular oscillations in the dominant high-frequency range (around 22 s), and slower oscillations (around 1 min) in the low-frequency range. Finally, expression of ACA2 with changes in four putative Ca2+ binding residues (ACA2{Delta}Ca2+) failed to rescue the RH growth phenotype in the aca2-2 mutant. Collectively, our findings indicate that ER-localized ACA2 and possibly ACA7 are crucial for modulating cytoplasmic Ca2+ signals, possibly composing a critical part of a negative feedback loop, and their absence leads to impairments in RH cell elongation.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

Polyploidization can profoundly affect plant metabolite biosynthesis, yet its influence on chemodiversity remains poorly understood, despite the central role of chemodiversity in mediating plant interactions with the environment. The coexistence of facultative apomictic and sexual reproductive systems across ploidy levels in Hypericum provides an excellent model for investigating the evolution of chemodiversity following polyploidization. We analyzed ploidy levels and leaf metabolic fingerprints across selected populations of three Hypericum taxa, H. maculatum, H. perforatum subsp. perforatum and H. perforatum subsp. veronense. Polyploidization was common across all three taxa. Leaf metabolic fingerprints were more pronouncedly differentiated by the ploidy level of the mother plant (F0) than that of the offspring (F1). Although unique metabolic features emerged in plants of most ploidy levels, diploid plants exhibited fewer metabolic features than polyploid plants. Higher Shannon diversity, functional Hill diversity, and intensities of features belonging to specific chemical families were associated with higher F0 ploidy levels in H. perforatum subsp. perforatum, but not in H. maculatum and H. perforatum subsp. veronense. Our findings demonstrate that polyploidization can lead to rapid shifts in chemodiversity across generations in Hypericum. The fast divergence in chemodiversity associated with polyploidization in H. perforatum may contribute to its remarkable invasive potential.

Photosynthesis is among the most consequential yet genetically complex traits in crop plants, and translating its natural variation into actionable genomic targets remains a central challenge for breeding climate-resilient varieties. To start addressing this, researchers are generating increasingly large, multi-environment field photosynthesis datasets. Yet, these data have been structurally under-analysed since their inception. Here we report the outcomes of the first dedicated hackathon focused on computational mining of such field data held in Accra, Ghana, in March 2026. Bringing together data scientists, plant physiologists, geneticists, and breeders from Europe and Africa, these interdisciplinary teams used photosynthetic data collected with hand-held fluorometers to genome-wide marker data across four crop species: cowpea (Vigna unguiculata), barley (Hordeum vulgare), common bean (Phaseolus vulgaris), and potato (Solanum tuberosum). Despite using different species and methods, independent teams identified the same three key findings. First, mechanism-informed feature engineering and dynamic modelling recover genetic signals that are not detected or discarded in standard analysis pipelines, resulting in traits with improved heritability and meaningful associations with yield. Secondly, machine learning methods proved effective at uncovering genetic associations, with temporally resolved features substantially outperforming single time-point measurements. Third, raw chlorophyll fluorescence and absorbance traces consistently contained more information and predictive power than the extracted parameters currently used. A defining feature of this event was having experimentalists and data scientists working together, enabling AI approaches to be grounded in domain knowledge and biological mechanisms rather than relying on data alone.

Tillering plasticity is a key adaptive trait in sorghum influencing resource use efficiency via a plants ability to adjust branching to neighbour density. Neighbour detection through red:far-red (R:FR) light sensing regulates this plasticity. While molecular pathways regulating tiller outgrowth are partly known, the genetic architecture underlying density-responsive tillering has not been resolved in any grass species. A sorghum diversity panel (n = 895) was evaluated over two growing seasons (2023 and 2024) with plant spacing ranging from 5 to 60 cm. A linear mixed model incorporating neighbour distance and tiller counts estimated genotype-specific response. GWAS was conducted on isolated plants (no neighbours within 60 cm) and on estimated responsiveness to neighbours. GWAS identified 52 baseline tillering QTLs and 50 for spacing responsiveness, with 10 overlapping, suggesting shared genetic control. Comparison with 41 R:FR pathway candidate genes revealed enrichment in responsiveness QTLs (5/50, 10%) versus baseline (0/52, 0%) (Fishers exact test, P = 0.025). Our model identified 40 unique density-responsive tillering QTL regions. Reducing genotype response to neighbour absence could be a selection target to develop water-efficient sorghum varieties where controlled architecture may be more valuable than natural plasticity.

Water-use efficiency (WUE), the ratio of accumulated plant biomass to water lost through transpiration has conventionally been determined using a destructive single-point measurement. Recent advances in high-throughput phenotyping now enable repeated, non-destructive estimation of biomass and WUE. However, these digital measurements must be statistically validated against conventional destructive methods to validate their use as reliable proxies. Therefore, we compared digital biomass determined point clouds produced from multispectral camera scanners with destructive harvests across eight harvests using Samsun tobacco grown under both drought and high-water conditions. WUE efficiency, calculated using the digital biomass estimated from a point cloud and gravimetric water use determinations, were compared to destructive harvest determinations. The coefficient of variation (CV) showed there were no significant differences in digital and destructive measurements for either biomass or WUE. Indicating that digital measurements can be used in place of destructive measurements. Drought plants used significantly less water and were significantly smaller than high-water plants from Harvests 4 through 8. However, there were no significant differences in the ratio of evapotranspiration to leaf area or WUE, indicating that drought plants were simply smaller and used less water than the high-water plants. This work validates that estimating plant biomass from a digital point coupled with continuous gravimetric determination of water use provides a reliable nondestructive measure of WUE in high-throughput measurements across the full plant life cycle.

Chloroplasts comprise photosynthesis and other important processes. Plants protect chloroplasts from stresses including Cd accumulation. Common terrestrial plants, excluders apply a set of mechanisms to restrict Cd penetration to chloroplasts. Removal of accumulated Cd from chloroplasts should also be a beneficial strategy. However, we do not know whether excluder plant species have ability to remove Cd from chloroplasts. We used barley as a common excluder plant species. To barley plants, we applied a model with two stable isotopes 111Cd and 114Cd to distinguish Cd accumulated earlier and later. A portion of Cd absorbed by roots continued translocation to shoot for some days after the external source of Cd was changed from one isotope to another. Chloroplasts acquired new portions of Cd and lost part of Cd accumulated earlier; a total Cd content remained rather unchanged. Cd loss from thylakoids was detected in vivo and in vitro. Cd loss from stroma and envelope was observed in vivo but not in vitro. Therefore, barley has at least two distinct mechanisms for Cd removal from chloroplasts: one from thylakoids and another from stroma. We hypothesized diverse chlorophagy pathways as a potential mechanism for Cd removal from chloroplasts. Cd accumulation by chloroplasts was mainly light-independent. In chloroplasts, Cd accumulated in vivo was tightly bound and mainly located in thylakoids. In vitro, chloroplasts from Cd-treated plants accumulated much less Cd than chloroplasts from untreated plants in a previous study. This implies reorganization of transport across chloroplast envelope membranes. HighlightsO_LICd was removed from thylakoids both in vivo and in vitro C_LIO_LICd was removed from stroma and envelope in vivo but not in vitro C_LIO_LIIn chloroplasts, Cd accumulated in vivo was tightly bound C_LIO_LICd accumulation by chloroplasts was mainly light-independent C_LIO_LIRoot barrier slowed down Cd translocation to shoot but not halted it C_LI

Root phenotypic plasticity is often proposed to improve crop performance under stress, yet it remains unclear how much plasticity is beneficial and whether adaptive responses require changes across many traits or adjustments in few specific traits. Using public data of 6,500 field-grown maize and barley plants, this study examined the extent and distribution of root plasticity, and when it is associated with yield stability. We quantified root plasticity across nine anatomical and architectural traits using complementary statistical models and applied a feature-discovery framework to identify the drought-associated optimal integrated phenotypes and determine whether plasticity toward these phenotypes improved yield stability. More plasticity did not mean greater yield stability. Neither the number of plastic traits nor the magnitude of plastic responses predicted yield stability. Rather, we identified species-specific high-yielding, stable integrated phenotypes defined by distinct trait configurations. Critically, genotypes whose plastic responses moved their root phenotype toward these targets achieved greater yield stability, whereas movement away from them was associated with lower stability. Root plasticity is adaptive when it shifts root phenotypes towards an optimal integrated phenotype. These findings show that the value of plasticity depends on the trajectory of phenotypic change rather than its magnitude alone.

Mechanistic models of plant regulatory networks typically require extensive parameterization, limiting their generalisation and scalability. Here we present a parameter-free, topology-driven model of shoot branching that predicts phenotypic outcomes from network structure alone. We constructed a signed, directed causal network by distilling regulatory relationships from the published literature spanning many laboratories, species, years, data types, and methodological frameworks. This extracted the essential logic of the system, consistent with developmental-biological reasoning and anchored in empirical evidence. Using PSoup, which automatically translates network topology into algebraic equations, the model propagates information across the network and predicts the qualitative direction of change relative to a defined baseline, mirroring the comparative framework of biological experiments. The pipeline, from network construction through automated equation generation to prediction, is transparent and reproducible. Trained against branching phenotype data with 78 diverse perturbations spanning genetic mutations and hormone treatments, the model achieved 86% accuracy in predicting branching direction. On an independent test set of 84 perturbations measuring bud release and gene expression at nodes not used during training, accuracy reached 75%. The approach highlighted deficiencies in our understanding of the topology of the network around SMXL 6/7/8 and ABA nodes. Other errors came mainly from modelling choices, such as the threshold for scoring a node as changed relative to baseline. Beyond shoot branching, this work demonstrates a general strategy for synthesizing biological knowledge into validated predictive networks, providing a foundation for both applied breeding and the advancement of fundamental biology.

1. Plants produce a great chemodiversity, which is the diversity of specialized metabolites (SMs). These SMs are produced in complex metabolic pathways and play an important role in inter-species interactions. There are numerous hypotheses about the evolutionary processes which brought about and maintain chemodiversity. Some have been partially tested in lab and field studies. However, some of their assumptions and predictions are better tested by quantitative modeling, and so far no quantitative model has investigated the role of metabolic pathways. 2. To close this gap, we developed an individual-based model for metabolic pathway evolution. It models enzymes creating metabolites with various modifications. Enzymes undergo inheritance and mutation. We used the model to compare the screening and interaction diversity hypotheses. 3. The screening hypothesis predicts promiscuous enzymes, genetic drift, the presence of many non-beneficial metabolites, and high metabolite richness. The interaction diversity hypothesis predicts specialized enzymes, selection, the almost exclusive presence of beneficial metabolites, and situation- dependent metabolite richness. We found that the patterns predicted by the screening hypothesis did not occur, while those predicted by the interaction diversity hypothesis did. 4. This provides reason to favor the interaction diversity hypothesis over the screening hypothesis when connecting empirical results to their evolutionary context

BackgroundThe capitulum of Asteraceae is a highly specialized inflorescence whose formation requires the coordinated regulation of multiple developmental processes, including floral organ identity and floral meristem determinacy. The LEAFY (LFY)-UNUSUAL FLORAL ORGANS (UFO) regulatory module is known to play an important role in flower development; however, naturally occurring mutations affecting this pathway have not been genetically characterized in gerbera (Gerbera hybrida). ResultsIn this study, we characterized a novel gerbera mutant identified during a commercial crossing program and named it marimo based on its green, spherical capitulum. Morphological observations revealed the repeated formation of secondary and tertiary floret-like organs within primary floret-like organs. Scanning electron microscopy showed that the epidermal structure of the green organs in marimo was similar to that of wild-type involucral bracts. RNA sequencing identified numerous differentially expressed genes between marimo and the wild type, and network and Gene Ontology analyses highlighted gene groups associated with flower development, reproductive organ differentiation, and tissue structure formation. RNA-seq analysis showed increased expression of LFY and reduced expression of GGLO1, a PISTILLATA/GLOBOSA-like B-class MADS-box gene, in the marimo mutant. RT-qPCR analysis of a segregating population further confirmed reduced GGLO1 expression in marimo-type individuals. In addition, a single-nucleotide deletion was identified in the coding region of UFO. This deletion was predicted to cause a frameshift and a premature stop codon. In selfed progeny of No. 251, the UFO genotype was fully associated with capitulum phenotype, and only individuals homozygous for the mutant allele exhibited the marimo phenotype. ConclusionsThese results indicate that the naturally occurring frameshift mutation in UFO is the strongest candidate variant underlying the marimo phenotype. RNA-seq analysis showed increased LFY expression and markedly reduced GGLO1 expression in the marimo mutant. Reduced activity of the LFY-UFO regulatory module may therefore have altered the expression of GGLO1 and other floral organ development-related genes despite the continued expression of LFY. These changes may have affected both floral organ identity and floral meristem determinacy, resulting in the formation of green involucral bract-like organs and the repeated production of floret-like organs. The marimo mutant provides a useful genetic resource for investigating capitulum development in Asteraceae and may also serve as breeding material for introducing novel ornamental traits into gerbera.

Stomata are the pores on plant surface, and these tiny pores are responsible for the flow of gas between plants and atmosphere. Currently, what effects of the broad and continuous increase in stomatal density achieved via genetic engineering on plant growth and development remain poorly understood. The 9 Arabidopsis transgenic lines with increased stomatal density were acquired through overexpressing FSTOMAGEN (the homologs of STOMAGEN, which are in Flaveria). The intermediate stomatal density (SD) lines exhibited increased trend in biomass. Compared with the lines with low SD, the biomass of Arabidopsis lines with intermediate SD (484 mm-2) significantly increased. There was a positive and significant correlation between biomass and relative water content. Across these transgenic lines, only during the earlier phase of growth, the leaf area exhibited a gradually increased trend as stomatal density increased, and there was both a significant linear relationship between SD and leaf growth rate and a strong linear relationship between SD and leaf area. In contrast, a clear relationship during the later phase wasnt observed. Under lower growth light intensity, there was an increased trend of biomass from other lines to the lines with intermediate SD, and the photosynthetic rate and stomatal conductance of the intermediate line were significantly increased. This study reveals plant-growth alterations that correspond to broad and near-continuous increases in stomatal density achieved via genetic engineering. Our study sheds light on the prerequisites for elevated stomatal density achieved via genetic engineering to promote plant growth.

Plant root cellular architecture and cell wall composition influence plant productivity, stress resilience, biotic interactions, and potentially soil carbon accumulation. This study establishes comprehensive compositional parameters for roots of a lowland switchgrass genotype, DVR3. Root traits were analyzed in 12.5 cm depth segments, from Zone 1 near the surface to Zone 4 down to 50 cm. Mean abundance ({micro}g/mg) for major cell wall components included cellulose 470 {+/-} 20, xylose 250 {+/-} 20, lignin 170 {+/-} 15, and total suberin 35 {+/-} 5. Composition and cellular anatomy varied with depth, in a partially coordinated manner. Cross sections showed extensive aerenchyma in mature root regions despite greater root mass density, corresponding to abundant lignin and cellulose. Deep roots were enriched for pectin-associated traits, including arabinogalactan II, homogalacturonan, and arabinose-associated linkages. Suberin content did not vary significantly, though Casparian strip formation, endoderm and exoderm thickening, and suberin surface staining progressed with development. Similar trends in root lignin and specific root length were observed for another lowland switchgrass genotype, AP13. These results suggest that it may be possible to genetically enhance native switchgrass root chemistry to promote soil penetration and below-ground carbon accumulation by reducing variability with development, potentially via cell-type specific adjustments. HighlightOlder, shallower switchgrass crown roots are enriched in lignin and cellulose, and deeper, younger roots are pectin-rich with juvenile cellular anatomy. A more uniform compositional distribution might enhance below-ground traits. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=120 SRC="FIGDIR/small/744798v1_ufig1.gif" ALT="Figure 1"> View larger version (47K): org.highwire.dtl.DTLVardef@159de98org.highwire.dtl.DTLVardef@124d714org.highwire.dtl.DTLVardef@1a49c14org.highwire.dtl.DTLVardef@2fa67_HPS_FORMAT_FIGEXP M_FIG C_FIG Schematic summary of switchgrass root anatomy and composition across four 12.5-cm depth zones of a 50-cm root system. Zone 1 represents older, shallow roots and Zone 4 includes younger roots and root tips. Representative cross-sections show greater aerenchyma development in older roots than in young root tips. The compositional heatmap shows higher cellulose, lignin, and xylose in Zone 1, higher pectin and nitrogen in Zone 4, and relatively little variation in suberin across zones.

Utricularia gibba L. is an aquatic carnivorous plant with a diverse set of capabilities. Reddening of traps frequently occurs in old in vitro cultures. While anthocyanins are often responsible for red coloration in plants, not every plant turns red. Stress factors like high light or excess sucrose have previously been shown to induce the formation of anthocyanins. Here, we hypothesized the red trap formation to be dependent on nutrient deprivation and tested nitrogen deprivation. The results suggest, that only in combination with light, nitrogen deficiency leads to the activation of the complete anthocyanin biosynthesis pathway and visible red coloration. However, in darkness, anthocyanin biosynthesis appears generally less active compared to light conditions and expression of most anthocyanin biosynthesis genes is not significantly upregulated under nitrogen deficiency.